Pulsar kicks via spin-1 color superconductivity 
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We propose a new neutrino propulsion mechanism for neutron stars which can lead to strong 
velocity kicks, needed to explain the observed bimodal velocity distribution of pulsars. The spatial 
asymmetry in the neutrino emission is naturally provided by a stellar core containing spin-1 color- 
superconducting quark matter in the A phase. The neutrino propulsion mechanism switches on when 
the stellar core temperature drops below the transition temperature of this phase. 
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Introduction. — The first pulsar was discovered more 
than 37 years ago 0. Since then, pulsars remain among 
the most interesting celestial objects in our Galaxy. Their 
observation allows to test experimentally the existence of 
the gravitational radiation as predicted by general rela- 
tivity Q ■ The recent discovery of a double-pulsar system 
Q is likely to further strengthen the status of pulsars 
as an astrophysical laboratory for testing general relativ- 
ity. In this Letter, we suggest that the physics of pulsars 
may also serve as the ultimate laboratory for testing the 
theory of strong interactions, quantum chromodynamics 
(QCD), in the regime of large baryon density. 

Matter at large baryon density is expected to be de- 
confined and color-superconducting (for reviews see, e.g., 
Ref. 0). Color-superconducting quark matter could exist 
inside neutron stars, whose central densities are highest 
in Nature. Therefore, it is important to study the physi- 
cal implications of such a possibility in detail. At present, 
this is not easy because our current knowledge regarding 
the ground state of neutral, /3-equilibrated dense matter 
is very limited [aJfilL Many color superconducting phases 
were proposed 0, M E3 EH EH E3- but it is not clear 
which of them can be realized inside stars. 

It has been known for a long time that typical spatial 
velocities of pulsars are an order of magnitude larger than 
those of their progenitors 0, . Taking into account 
the violent conditions at the pulsar birth, this may not be 
so surprising. Even a small asymmetry in the supernova 
explosion may result in a kick velocity of several hundred 
km s^ 1 ll5J|. However, the bimodal velocity distribution 
of pulsars |lfl ll^ , which is unlikely to result from a sin- 
gle physical mechanism, might be more surprising. If one 
associates the low- velocity component (< 100 km s _1 ) 
with asymmetric supernova explosions, then what is the 
origin of the high-velocity component (> 500 km s -1 )? 
Several mechanisms were proposed jig, llSl l2lj , but the 
issue does not seem to be settled. For short reviews see, 
for example, Ref. |2l| . 

In this Letter, wc propose a new neutrino propulsion 
mechanism, resulting from a color superconductor in the 
transverse A phase |22j. In this phase, quarks of the same 
flavor form Cooper pairs with total spin one. The neu- 



trino emission from this phase, dominated by direct Urca- 
type processes, is not symmetric in space. This emission, 
as we shall see, provides a natural mechanism to power 
strong, e.g., of order 1000 km s _1 , velocity kicks for neu- 
tron stars, which could explain the high-velocity com- 
ponent in the pulsar distribution 0, • Unlike most 
other mechanisms, this one turns on not immediately af- 
ter the supernova explosion, but after the temperature of 
the stellar core drops below the critical temperature of 
the A phase. A distinctive prediction of this mechanism 
is the alignment of the kick velocity direction with the 
rotational axis. 

Neutrino emission. — Let us start by outlining the main 
steps in the derivation of a general expression for the neu- 
trino emissivity in spin-1 color superconducting phases. 
We use the Kadanoff-Baym formalism El to de- 
rive the following differential expression for the emissivity 
Hi: 
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where Gp is the Fermi coupling constant, fi e is the elec- 
tron chemical potential, and 5P* = P* — <5g^ e . Here, 
particle four-momenta are denoted by capital Latin let- 
ters, while the absolute values of the three-momenta are 
denoted by lowercase letters. The metric tensor is g\ a = 
diag(l, —1,-1,-1). The Bose and the Fermi distribution 
functions are denoted by = [exp(w/T) — 1] _1 and 

uf(oj) = [exp(w/T) + respectively. The lepton ten- 
sor L\ a (P e ,Pv) is defined as follows: 

Lx4Pe, P.) = Tr [P e K 7 «7a(l - 7 5 R P 7p7a(1 - 7 5 )] ■ 

(2) 

Finally, the last factor in the integrand on the right hand 
side of Eq. is the imaginary part of the retarded po- 
larization tensor of the W-boson, 

n ACT (Q) = 2 £ / 72^3 Tr [r x S(K)ris(K + q)] , 

(3) 
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with the trace running over flavor, color, Dirac and 
Nambu-Gorkov indices. The quark propagator S(K) is 
diagonal in flavor space, and its components have the 
following Nambu-Gorkov structure: 
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We consider the ultrarelativistic limit. The explicit color 
and Dirac structure of and 2^ can be found in 
Ref. |22. Here, we note only that the poles of the quark 
propagators appear at fc = k + (j,f (antiquarks) and at 



rj 



^{k-tijf + Xk^fW r=l,2,3, (5) 



where </>/ is the gap parameter, and the functions Xk^r 
are specified by the choice of the phase. 

The order parameter of the A phase has a special direc- 
tion in color space: quarks of one color do not pair. Also, 
it has a special direction in momentum space, say, the 
z-direction. If 6-^ denotes the angle between the three- 
momentum of a quasiparticle and the z-axis, the three 
low-energy quasiparticle modes in the transverse A phase 
are defined by Afe,i = (1+| cos#k|) 2 , A/-^ = (1 — I cos6>k|) 2 , 
and Afe.3 = 0. Here, "transverse" refers to the fact that 
quarks of opposite chirality form Cooper pairs. 

The explicit expressions of the vertices in Eq. J3J read 
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f(v,o= E 

ei,e 2 =± q 



1 



dx u dx d 



and 



( e -ei^ + 1) ( e e 2 £~ d + 1) ' 

(9) 

u, d. Note that 



e f = ^jxj + (1 + with / 

< G 0„, ip d ) < 1 and G(0, 0) = 1. 

In passing, we note that the pair breaking processes 
[i^ do not play any significant role in the case of spin-1 
color-superconducting quark matter under consideration. 
The corresponding contribution to the emissivity |29| is 
parametrically suppressed by factor T/fi e ~ 10~ 3 . 

From Eq. JJJ, we see that the (anti-)neutrino emissivity 
in the A phase differs only by a factor of order 1 from the 
corresponding result in the normal phase of quark matter 
|26|,|27j. Therefore, the A phase should have qualitatively 
the same effect on cooling of stars as the normal phase. 

Nevertheless, the neutrino emission from the A phase 
is very unusual. It is not symmetric with respect to re- 
versing the direction of the z-axis. To quantify the asym- 
metry, we calculate the value of the z-component of the 
momentum carried away by neutrinos per unit volume 
of quark matter, per unit time. This is obtained by re- 
placing one power of p v on the right hand side of Eq. Q 
by p u cos 8 Pv . Taking into account that the neutrino and 
the antincutrino emissions give the same contributions, 
we arrive at the final result 251 
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where the flavor matrix t± = (ti ± ir^jl is constructed 
from Pauli matrices. 

By substituting the quark propagators J3J and the ver- 
tices © into Eq. ©, we calculate the imaginary part of 
the retarded polarization tensor. Then, by making use of 
the result, we derive an expression for the emissivity in 
the following approximate form [25j : 
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where e£ r ^ is the partial contribution involving the rth 
type quasiparticle modes, see Eq. (0. The factor 2 comes 
from taking into account both the neutrino and the anti- 
neutrino emissivities. In the final result, the contribution 
of the ungapped modes v v is, up to a factor 1/3, the 
same as in the normal phase of quark matter |2(| [2?| ■ 
The contribution of the other two modes is suppressed 
by the following function: 



G(ip u ,ip d ) 
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where 
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Note that H(0, 0) = which is consistent with the fact 
that the momentum kick is vanishing in the normal phase 
of quark matter. The numerical result for the function H 
at equal values of its two arguments is well approximated 
by the following expression: 
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with hi = 0.3068, h 2 



-0.1977, h 3 = -0.7838, h 4 



1.0286, h 5 = -0.3539, and r = e^/n w 0.8125 is 
the ratio of the critical temperature to the value of the 
gap at T = in the A phase, i.e., tq — T c /4>$ which 
is expressed in terms of the Euler constant 7 ~ 0.577 
and C = In 2 - 1/3 0. Note that J2l=i h n = 0. The 
representation in Eq. (|12fl is particularly convenient when 
the following simplified temperature dependence for the 
gap parameter is used: 4>(T) = <j>Qy/l — (T/T c ) 2 . In this 
case, the function H(T) becomes a polynomial: H(T) = 
En=l K{T/T c ) n for T < T e , and H{T) = for T > T c . 
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The physical reason for the breakdown of the reflection 
symmetry lies in the pairing pattern of the transverse A 
phase. Quasiquarks of the first branch, r — 1, have helic- 
ity +1 if the projection of their momentum onto the z- 
axis is negative, cos#k < 0, and helicity —1 if cos6>k > 0. 
Quasiquarks of the second branch, r = 2, have oppo- 
site helicities. Only left-handed quarks (in the ultrarel- 
ativistic limit, quarks with negative helicity) participate 
in the Urea processes. Thus, the quasiquarks of the first 
(second) branch contribute only if their momenta are in 
the upper (lower) half-space. Taking this into account, 
the effective branch relevant for the emission has gap 
4>cS ~ l + cos#k, which discriminates between +z and — z 
directions, see Fig.Q] Since neutrinos are emitted prefer- 
ably in the direction opposite to the quark momenta, this 
asymmetry manifests itself in the neutrino emission. 

Calculation of the velocity kick. — In order to make a 
simple estimate of the velocity kick due to the neutrino 
emission from a color-superconducting quark matter core 
in the A phase, we use the following model time de- 
pendence of the core temperature (see, e.g., Ref. [27^: 
T(t) = T (t /t)^ 4 . This has the power dependence which 
is characteristic for bulk matter cooling by neutrino emis- 
sion with e„ ~ T 6 , see Eq. Q, provided the specific heat 
of matter is cy ~ T. Here, we assume that T c < To, i.e., 
the system is too hot for spin-1 pairing initially, and then 
it cools through the transition point. 

The velocity kick for a star of mass IAMq with the 
quark core of radius R c is given by the expression: 
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where we used the approximate expression in Eq. (|12H . 
as well as a simplified temperature dependence of the 
gap parameters, </> u (T) = ^ d {T) = 4> ^1 - (T/T c ) 2 . The 
notation t c = to(To/T c ) 4 stands for the time when the 
temperature of the quark core drops below the critical 
value. 




From Eq. (|13|l . we can also derive the expression for 
the maximum velocity kick (t — oo): 
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The results for the maximum velocity kicks are pre- 
sented graphically in Fig. [21 To make the plot, we take 
a s = 1, and use the initial condition: To = 100 keV at 
to = 100 yr. Also, we use the following values of the 
chemical potentials: /i n = 400 MeV, Hd = 500 MeV, 
fi e = 100 MeV. 

When the stellar age t is finite, all curves in Fig.|3shift 
to the right. Numerically, the value of the shift is 10 keV 
at t = 10 6 yr and about 30 keV at t = 10 4 yr. This is easy 
to understand: at t < oo, the step function in Eq. i|13|) 
cuts out a range of low values of critical temperatures, 
T c < To(to/t) 1/>4 , which arc not yet accessible at time t. 
Other than this shift, the shape of the curves remains 
similar even at t as low as 10 3 yr. 

Discussion. — The main observation of this Letter 
is that the neutrino emission from the spin-1 color- 
superconducting transverse A phase (unlike from other 
spin-1 color-superconducting phases 25]) is not symmet- 
ric with respect to reversing one spatial direction. 

Therefore, we propose a hypothesis according to which 
the inner cores of some (hybrid) neutron stars are made 
of the A phase. The neutrino emission from such stars 
could generate strong velocity kicks, needed to explain 
the observed bimodal velocity distribution of pulsars 
[lil Il7| . Its bimodal structure results from an overlap 
of two distributions: one describing "normal" neutron 
stars, and the other describing hybrid stars with color- 
superconducting quark cores in the A phase. 

The above hypothesis has several specific (and, thus, 
falsifiable) predictions, which are directly related to the 
nature of the A phase. Let us start by discussing the 
direction of the kick velocity. As should be clear from our 
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FIG. 1: Gap functions for the first (left) and the second (mid- 
dle) excitation branch with specified helicities of quasiparti- 
cles in the upper and the lower half-spaces. The "effective" 
gap relevant for the neutrino emission is shown on the right. 



FIG. 2: The dependence of the maximum velocity kick of a 
neutron star versus the value of the critical temperature of 
the spin-1 color superconducting phase transition. Results for 
several values of the quark core radius are shown. 
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analysis, this is aligned with the direction spontaneously 
picked by the spin-1 condensate. In the case of an infinite 
isotropic medium, this direction can be arbitrary. Inside 
a rotating and magnetized star, however, the degeneracy 
with respect to the orientation is likely to be removed 
because of a nonvanishing interaction with the angular 
momentum and/or with the magnetic field. Therefore, 
the kick velocity averaged over time should be parallel 
to the rotational axis. The observational data seems to 
favor this possibility [sol ]. 

The proposed neutrino propulsion mechanism also pre- 
dicts a correlation between the velocity distribution and 
the age of young neutron stars. This has not been seen 
in the observational data yet ^(|. Most likely, however, 
this correlation is hard to detect because it takes typi- 
cally less than 10 4 yr for a star to accumulate almost all 
of the maximum velocity kick. The future high statistics 
studies of young neutron stars may resolve the issue. 

Along with the advantages of the proposed mechanism, 
it might be also appropriate to mention several potential 
difficulties. One of them is a rigorous justification that 
the spin-1 color-superconducting A phase is the ground 
state of dense quark matter inside a star. In the case 
of an infinite isotropic medium, it has been shown that 
the color-spin locked phase is the favored spin-1 phase at 
asymptotically large density 0,H3- We could only spec- 
ulate that this may change when the density is realistic, 
and the rotation and the magnetic field are accounted 
for. The other potential difficulty is related to fast cool- 
ing of the A phase by direct Urea processes, which might 
not be compatible with the observational data [sil ]. 

The mechanism, proposed here, may well have other 
important implications that we did not discuss in this 
Letter. We hope, however, that they will be addressed in 
the future studies, providing much stronger tests for the 
hypothesis of the hybrid stars accelerated by the asym- 
metric neutrino emission from the A phase. If it passes 
the tests, the observational data from pulsars could shed 
light on some details of the QCD phase diagram. Thus, 
the pulsars should be viewed not only as a unique labora- 
tory for testing the theory of general relativity, but also 
the ultimate laboratory for testing the theory of QCD. 
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The model time dependence of the temperature, 

T(t) 



To 



1/4 



(1) 



which was used in our calculation, is unphysical. As a result, the estimate of the kick velocities due to spin-one color 
superconductivity was wrong. 

The correct time dependence is approximately given by 



T(t) = To 
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which is the same as for normal quark matter. This time dependence is obtained by integrating the energy balance 
equation, e„ = —cydT/dt, in which an approximate analytical expression for the specific heat of quark matter, 
cy = (/i 2 + (J%)T (sec, for example, Ref. Q) and Iwamoto's analytical expression for the neutrino emissivity e v Q 
were used. This is justified for the long-term stellar cooling after the thermal relaxation epoch is completed (t > 10 2 yr) 
and before the surface cooling by photons starts to dominate (t < 10 5 yr), e.g., see Ref. Q. In the above expression, 
t is defined by 
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Taking this into account, Eqs. (13) and (14) in our Letter should be replaced by 



Sv 



(tot) 



I AM, 



© 



457a, 
945 



and 



G p^eHu^d 



Sv„ 



4tt R 3 c 
T 1.4M 



^0 



-T£T 2 c T6{t-t c ) 



5 

71=1 



t-to+T 



(2+n)/4' 



0.033 a s G%fj, e p,ufJ-d 



4tt R 3 c 
YlAM, 



1 1 c T i 







(3) 



(4) 



(5) 



respectively. Here t c = t(Tq/T* — 1) + to is the time when the temperature is equal to the critical value. 

The revised expression for 6v max differs by a factor r/to w 10~ 7 from the old one. Consequently, the predicted kick 
velocities become extremely small. For example, instead of velocities of the order of <5i> max ~ 1000 km/s, one gets 

fcmax ~ 10~ 4 km/s. 
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